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Abstract 

By use of the multi-Josephson junction model, we investigate voltage-biased I-V 
characteristics. Differently from the case of the single junction, I-V characteristics 
show a complicated behavior due to inter-layer couplings among superconducting 
phase differences mediated by the charging effect. We show that there exist three 
characteristic regions, which are identified by jumps and cusps in the I-V curve. In 
the low voltage region, the total current is periodic with trigonometric functional 
increases and rapid drops. Then a kind of chaotic region is followed. Above certain 
voltage, the total current behaves with a simple harmonic oscillation and the I-V 
characteristics form a multi-branch structure as in the current-biased case. The 
above behavior is the result of the inter-layer coupling, and may be used to confirm 
the inter-layer coupling mechanism of the formation of hysteresis branches. 

1 Introduction 

I-V characteristics of the high T^ superconductor Bi2Sr2CaiC208 (BSCCO) shows a strong 
hysteresis, producing multi-branches |, ||. Many experiments have been reported to 
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investigate properties of hysteresis branches ^ ^, ^, ^ |I^, and it has turned out 
that there is rich physics in properties of I-V characteristics. Substructures are found in 
hysteresis branches, and are discussed in relation with phonon effects The d-wave 
effect is also discussed in the higher voltage region, where the tunneling current becomes 
important [|, p^ . 

It is well known in a single Josephson junction that the application of a constant 
voltage-bias leads a rather simple phenomenon of the harmonic current-oscillation with 
the frequency proportional to the applied voltage This is due to the fact that 

the time-derivative of the phase difference is directly connected with the voltage, and 
that there is no room for the phase to move freely. Therefore, the most of the above 
mentioned experiments were performed with the current biased arrangement, aiming to 
investigate properties of I-V hysteresis branches, and, as long as we know, no systematic 
experiment with the voltage bias has been reported. We will point out in this paper 
that anomalous behaviors are expected, theoretically, in the low voltage-biased region in 
the case of the multi-layered intrinsic Josephson junctions of high Tc superconductors, 
and that the voltage-biased experiment gives additional information on behaviors of I-V 
characteristics. 

In the multi-layered intrinsic Josephson junctions of high T,, superconductors, inter- 
layer couplings among superconducting phase differences become important due to the 
atomic scale of the thickness of the superconducting layer. In the ref. [jl2| the mechanism 
of the inter-layer coupling is proposed as the charging effect of superconducting layers; the 
charging effect induces the variation of electric field and mediates the inter-layer coupling 
among superconducting phase differences. It has been shown that, by this mechanism, a 
branch structure in the I-V characteristics is obtained as the intrinsic nature of the strong 
anisotropy of high T^ superconductors ||T^. In the previous paper we have investi- 
gated the origin of hysteresis branches, and have shown that hysteresis jumps are induced 
when solutions change in nonlinear-coupled equations of phase-differences. Namely, the 
number of rotating phases and its patterns classifies solutions. At a point of a hysteresis 
jump, its number or distribution pattern is changed. 

The main difference between the cases of the single junction and the multi-junction lies 
in the fact that the applied voltage gives restriction only on the total voltage. Therefore, 
each junction may behave in a self-adjusted way, if an interaction effect allows it. In a 
simple array of single junctions, each junction feels only the A^-th of the applied voltage, 
and one cannot expect much difference from the case of the single junction. However, in 
the intrinsic Josephson array of the high T^ superconductors, the existence of an inter- 
layer coupling affects the collective motion of phase differences. Furthermore , the voltage 
control and current control lead different situations. Specially, in the voltage control case, 
one can reaches transient unstable states of the current biased case, since the applied 
voltage still enforces the total behavior of phase differences. In this paper, we investigate 
I-V characteristics in voltage-biased cases of high T^ superconductors. Since there is no 
experiment available to refer at present, several cases with possible surface conditions will 
be presented. 
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In the next section, we summarize the formula for the voltage-biased case in the 
multi-Josephson junction model. In Sect. 3, we present numerically simulated results, 
and present I-V characteristics of the voltage-biased case. It will be shown that certain 
anomalous behavior of I-V characteristics appears in the applied low voltage region. We 
show that the voltage region is divided roughly into three regions, an anomalous periodic 
region, a chaotic region and a region of hysteresis branches. In an anomalous periodic 
region, the current shows a trigonometric functional increase and rapid drop. In a chaotic 
region, the total current behaves chaotically. In the last region, the behavior of the current 
is quasi-harmonic. Physical origin of these regions will be discussed. Sect. 4 is devoted 
to concluding remarks. 



Formula for the voltage biased multi-Josephson junc- 
tion 



In this section, we rewrite the formula for the multi-Josephson junction of ref. suitable 



to the voltage-biased analysis. The notations used here are same as those in ref. |14 . 
Namely, we consider the + 1 superconducting layers, numbered from to A^. We denote 
the gauge invariant phase difference of the (/ — l)-th and l-th superconducting layer by 
(p{l), and voltage by V{1). The widths of the insulating and superconducting layers are 
denoted by D and s, respectively. At the edges, the effective width of the superconducting 
layer may be extended due to the proximity effect into attached lead metals. The widths 
of the 0-th and A^-th superconducting layers are denoted by sq and sn, respectively. 
In the voltage-biased case, the total voltage is equal to the applied external voltage 

N 

Ve,t = T.Vl- (1) 
1=1 

In order to perform a numerical simulation under the condition (|l]), we rewrite the equa- 



tion for the total current J (Eq. (1) in ref. |]T4|) and the equation relating the time- 



derivative of phase- difference and voltage (Eq. (3) in ref. [|T4|) in the following way. 



First, we note that the total current J is expressed by the external voltage Vext as 



J 1 



Vext _|_ 1 ^ Vext j 



The current J is normalized by the critical current Jc, jc{l) = Jc{l)/Jc with Jc(0 being 
the critical current for the l-th junction. The time t is normalized by the inverse of the 



plasma frequency, ojp = y x^^T^' ^i^^ e being dielectric constant of the insulating layer. 
The voltage V{1) is normalized by \^ = The parameter /3 is given by /5 = We 
define the voltage difference v{l) as 

Vil) = ^Vext + v{l) . (3) 
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Then we have the following coupled differential equations, 



and 



1 dv{l) 
Up dt Vp 



v{l)_ 

Vr, 



N 



i'=i 



(4) 



1 ^ 




1 V,^t ^ v{l') 



N K 



(5) 



A dissipation effect [|1J, [T5|, 0in superconducting layers will be neglected in this paper, 
and the matrix A is given in Eq. (4) in ref. ||14|. It should be noted that Eq. (^) gives 



I Q N N 

-7^E-(0 = -/?E-(0. 

u}p at 



Therefore if 



N 

1=1 







(6) 



(7) 



is satisfied at the initial time, it is satisfied in all time. In the next section, we present 
results of solving Eqs. and (^) by numerical simulation. The restriction from the total 
voltage is achieved by putting the initial condition Eq. (0). The equation for v'(/), Eq. 
(I), indicates that the phase differences move according to -^Vext as average. However, 
due to the mutual interaction, each ip{l) shows a complicated behavior which affects the 
I-V characteristics and the behavior of the total current. 



3 Results of numerical simulation 

We have solved the coupled differential equations @) and (|^) by use of the fourth order 
Runge-Kutta method. The average current J is obtained by 

J = -—- / dtJ{t) . (8) 

-'-max min ''T^i^ 

We choose the parameters as = 10, a = 1.0 and (3 = 0.2, which has been used in ref. 
In the simulation we have chosen the time step Updt = 1.0 x 10^^, [3ujpTmin = 10.0, 



(3ujpTmax = 310.0 and the voltage step dVext = 0.01. 

As Eq. (^) shows, the phase increases according to / dtj^Vext when Vext is small. 
When (p{l) reaches the value |, the phase goes into a rotating mode as can be seen from 
the analogy of the motion of a pendulum with a sin(y9(/))-nonlinear term. This indicates 
that the whole motions of phases are much affected by the phase with the lowest jc- 

It may be difficult to arrange two electrode at edges exactly in the same condition. 
However, in order to show how surface conditions affect the behavior of I-V characteristics, 
we present the following three cases; l)Jc(0 symmetric but widths of surface are 
asymmetric, 2) both 7^(0 widths of surface sjTiimetric, and 3) Jc(0 asymmetric and 
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widths of surface symmetric. The case with Jc(/) asymmetric and widths of surface 
asymmetric shows similar behaviors as the case 3). It should be noted that, even if the 
asymmetry between sq and is introduced, the observed critical current is Jc in the 
current-biased case. 

We first assume that all jc{l) are equal. 



and choose, for example, an asymmetric boundary condition sq/s = 1.0 and sn/s — 2.0. 
In the present analysis, it is enough to have slight difference between sq and sjv, the 
difference of proximity effect. 

In Fig. 1(a), we show the overall behavior of I-V characteristics obtained from the 
adiabatic increase and decrease of the applied voltage (circle points) and those obtained 
from abrupt application of the voltage ( x points) . We can see that the I-V characteristics 
show many hysteresis jumps, forming a branch structure. In the abrupt application of 
voltage, the data points in the low voltage region are scattered, although they lie on some 
hnear I-V branches. Each I-V branch is classified by a number and pattern of rotating 
phases. When the voltage is large enough, all junctions have the rotating phase and 
the single I-V characteristic is obtained. In Fig. 1(a), we notice that there exists the 
region where the adiabatic and abrupt applications of voltage give the same results, in 
the low voltage region before the branch structure starts. In Fig. 1(b), we expand the low 
voltage region. We can identify the obvious branches labeled as Ii, I2, II, IIIi and III2. 
The crosses (x) are for the abrupt voltage application. The fines I and II are same for 
three cases of the adiabatic increase, decrease and abrupt appfication of the voltage. The 
branches IIIi and III2 are finear. The calculation shows that the I-V relations for IIIi and 
III2 are same as those obtained from the hysteresis branches in the I-V characteristics 
of the current-biased case. On the line IIIi, the phase (p{l) is rotating and on the line 
III2 the phases (/?(1) and V'(A^) are rotating. Other phases (/?(/) are oscillating. The cross 
points are grouped as one rotating-phase and two-rotating phase. The scattered points 
are identified by the distribution of the rotating-phases among N junctions. 

In Fig. 2, we show the I-V characteristics, when we choose a symmetric boundary 
condition sq/s = 1.0 and sn/s — 1.0. The branch II in Fig.l is separated into IIi and II2. 
The branch III2 is finear and is exactly same line as III2 in Fig. 1. Although the branch 
II2 is linear as average and lies on the same line as IIIi, the calculated points are largely 
scattered, indicating a certain kind of instability. In order to obtain Fig. 2, we have to 
ensure the symmetry. 



in each time step by taking the average of (/?(/) and Lp{N + 1 — /), for example. Otherwise, 
a very small round error rapidly develops in the region II, and we obtained the similar 
result of Fig. 1(b), although we start with the symmetric boundary condition. 

In order to see what kinds of solutions are realized in the regions I and II, we plot, in 
Fig. 3, the time-dependence of the phase (p{l) and the current J for the case of Fig. 2. Fig. 



jc{l) = 1.0 . 



(9) 



ip{l) ^ip{N + l-l) , V{1) = V{N + 1-1) , 



(10) 
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3(a),(b) are for the branches Ii , and Fig. 3(c),(d) are for IIi. Fig. 3(b) shows periodic 
current behaviors, while Fig. 3(d) shows chaotic current behaviors. The results show 
the following. All (p{l) increase slowly in the same way according to the applied voltage. 
When {p{l) and (p{N) at the edges reach the value |, they increase rapidly in the phase- 
rotating mode. Other </?(/) 's swing back. The phases ^p{l) and ip{N) damp their motion 
by the effect of resistance /?, and move slowly again in the same way as other phases. 
As the result, the current behaves periodically with trigonometric functional increase and 
rapid drop. The regions Ii and I2 are different in the values of phase rotation, 27r and 
4:71, before damping. In Fig. 3(c), not only the phases at the edges but also other phases 
(/?(/) go into the rotating mode. Since oscillation and rotation of phases are interface in 
a complicated way, the chaotic current behavior is produced. The regions IIi and II2 are 
different in the values of phase rotation, 27r and An, at the edges. 

The above results lead us the following understanding of solutions of the multi- 
Josephson junction in the voltage-biased case. The applied voltage Vext forces the phase- 
differences ip{l) of junctions to rotate. When some ip{l) reaches |, it goes into the natural 
rotating mode with other phases swinging back. How much it rotates depends on its speed 
and resistance, forming new non-ohmic branches. When the rotation and oscillation of 
phase differences start to interact, solutions go into a chaotic region. When the voltage 
becomes large enough for certain phase difference to keep rotating, the I-V characteristics 
become ohmic and form multi-branches. 

From the above understanding, we can immediately see that the I-V characteristics 
may change when there exists a certain junction with weaker jc- Most likely, the surface 
layers may be damaged in the preparation process of electrodes. In fact in the previous 
analysis |]14|, we have to choose reduced Jc values at the edges in order to fit the data. It 
was also reported that for thin film mesa structures the preparation process for the upper 
electrode can produce a cut of the first junctions, leading to an area about half size and 
thus half critical current for those junctions As an example, we choose 

Ul) = 0.5 , Ul) = 1.0 (/ ^ 1) . (11) 

The result is presented in Fig. 4. The chaotic region II is missing. In the inset Fig. 
4(b), we show the enlarged figure of the region I. There are some drops and quasi-periodic 
structure in the I-V characteristics. Each part has a different value of (y9(/)-phase rotation. 



27m. We note that, in ref. |T8[, the appearance of certain periodicity in I-V characteristics 
has been reported, though it is for the sample of submicron junctions with no branch 
structures. It is a future problem to see if there is some relation with the present analysis. 

We also point out that, in the analysis of Fig. 1, the regions I-II are obtained also 
when the external voltage is applied abruptly. Adiabatic and abrupt voltage application 
differ in the initial total energy. Therefore, the present result shows that the regions I-II 
have the single stationary state. 
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4 Conclusion 



In this paper, we have investigated the I-V characteristics of the voltage-biased case. 
Differently from the case of the single junction, I-V characteristics shows the characteristic 
behavior originated from the mutual inter-layer coupling mediated by the charging effect. 
Although each junction feels as average the external voltage as ^Vext, there is a still 
room for each phase to move self-consistently according to the interaction. Therefore, 
when the voltage is low, there is a competition of the voltage forced motion and the 
motion induced by the interactions. Adiabatic increase of the external voltage induces 
the following successive transition. In the low voltage, one has the region with a periodic 
oscillatory current of an asymmetric waveform, and then a possible region with a chaotic 
current oscillation follows. In these regions, non-ohmic branches are formed. In higher 
voltage, the region with a stable harmonic oscillatory current appears, having ohmic I-V 
branches. When there is a junction with a weaker jc, the associated phase difference is 
easier to go into the rotating mode. Then the I-V characteristics have a quasi-periodic 
behavior in its non-ohmic region. The present analysis shows that voltage biased I-V 
characteristics reveals the importance of inter-layer coupling more strongly, and gives 
additional information to the current biased case. 
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Fig. 1. I-V characteristics with an asymmetric boundary condition 
(a) for the wide range of V and (b) for the expanded low voltage region. The 
circles are for the adiabatic voltage change. The crosses are for the abrupt 
voltage application. Parameters are a = 1.0, (3 = 0.2, = 10, — = 1.0 and 
^ = 2.0. 
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Fig. 2. I-V characteristics witii a symmetric boundary condition 
Parameters are a = 1.0, (5 = 0.2, AT = 10, ^ = 1.0 and ^ = 1.0. 
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Fig. 3. Time dependence of (/?(/) and J/Jc for the branches Ii ( (a)-(b) ) and 
III ( (c)-(d) ) 

Parameters are a = 1.0, 3 = 0.2, = 10, ^ = 1.0 and ^ = 1.0. 
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Fig. 4. I-V characteristics with a jc(Z)-asymmetric boundary condition 
Parameters are a = 1.0, /3 = 0.2, N ^ 10, ^ ^ 1.0, ^ = 1.0, jc{0) = 0.5 and 
je(0 = l-0(/7^0). 
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